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explained by a linear response analysis. Using direct simulations and the perturbative 
analysis it is shown that in the presence of a two-colour supcrlattice the excitation 
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the quasi Bose-glass is the appearance of low-lying resonances and a suppression of the 
dominant resonance of the Mott-insulator phase. 
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1. Introduction 

Ultracold bosonic gases in optical lattices offer a supreme laboratory for the study 
of a wide range of phenomena in strongly correlated quantum systems. Due to the 
flexibility of the experimental methods available for preparing and controlling the 
relevant parameters, one can address fundamental questions of importance in condensed 
matter physics and quantum information processing. 

In addition to the powerful experimental tools for preparation and control, several 
techniques for probing the quantum state of these many-body systems are available. The 
simplest experimental observable is the matter-wave interference pattern after realease 
of the atoms from the lattice, which led to the observation of the superfiuid (SF) to 
Mott-insulator phase transition (MI) ^ 13 El IH El • Another tool for characterizing the 
different quantum phases are dynamical excitations. The excitation via temporal lattice 
modulations corresponding to Bragg spectroscopy O IZl provides an excellent probe 
for the observation of quantum phase transitions. In contrast to static perturbations 
generated, e.g., by tilting the optical lattice [HIHl, Bragg spectroscopy provides a very 
precise probe for the response at a well-defined excitation energy |3l l9l ITT] IT2j. 

Recent experiments also explore the effects of irregular lattice potentials on 
the phase diagram. One can employ speckle patterns to create random lattice 
potentials (THl (HI lEl HEl (HI or use a superposition of two standing-wave lattices of 
different wavelengths to generate a two-colour superlattice potential [IHl 112] • Recently, 
the response of the Bose gases in time-dependent superlattices was investigated 
experimentally by Fallani et al. ^7]. The spatial structure of the lattice potential gives 
rise to additional quantum phases, such as localized [201 ^^"^ (quasi) Bose-glass phases 
• The systematic experimental investigation of the resulting phase diagram requires 
powerful probes that allow to distinguish the different phases. In addition to the matter 
wave interference pattern, the dynamical behaviour of the system provides important 
information. 

The aim of this paper is to study the response of zero-temperature Bose gases 
in one and two-colour superlattices and to identify possible experimental hallmarks of 
the individual quantum phases. To this end, we simulate the time-evolution of the 
system in the presence of a time-dependent lattice potential in the framework of the 
Bose-Hubbard model. In section |2] we formulate the time-dependent Bose-Hubbard 
Hamiltonian and discuss the numerical approach used for the time-evolution including 
a physically motivated basis truncation scheme. In preparation of the investigation of 
superlattice potentials, we study the dynamic response of Bose gases in regular lattice 
potentials in section El In addition to the full time-dependent simulations, we use a 
linear perturbation analysis pTl IT^ to characterize and interpret the response of the 
system. In section jH we extend this discussion to two-colour superlattices and identify 
possible dynamical signatures of the Mott-insulator to quasi Bose-glass transition. 
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2. Bose-Hubbard model and time evolution 
2.1. Bose-Hubbard model 

The Bose-Hubbard model has proven to be the appropriate framework for the description 
of zero-temperature Bose gases in optical lattices jlj . It describes the full phase diagram 
ranging from weakly interacting superfiuid phases to the regime of strong correlations, 
e.g. in the Mott-insulator phase. Assuming that the lattice potential is sufficiently 
deep, the single-particle Hilbert space can be restricted to the lowest energy band. A 
suitable single-particle basis is given by the localized Wannier functions of the lowest 
band. A state of bosons at / lattice sites can be represented by an /-tuple of 
occupation numbers {ni, 77,2, ■ ■ ■ , n/} of the individual sites [110 CH]- The Fock states 
\{ni, 71.2, ■ ■ ■ , ni}a) for all possible compositions of occupation numbers span a basis of 
the single-band Hilbert space. 

The Bose-Hubbard Hamiltonian for a single-component Bose gas in one dimension 
reads jl]: 

H = - J ^ (aj+ia^ + h.a.) + ^ e^hi + — ^ n^ (iii - 1) , (1) 

i=l i=l ^ i=\ 

with creation (annihilation) operators a| {k^ ) for a boson at site i and occupation number 
operators The first term in accounts for the hopping between adjacent 

sites with the tunnelling matrix element J. The second term introduces a site-dependent 
single-particle energy which is used to describe, e.g., an external harmonic trapping 
potential or a superlattice potential IT^ . The third term of the Hamiltonian (Q) 
accounts for the on-site two-body interaction of the atoms with interaction strength V . 
Throughout this paper we use periodic boundary conditions, i.e., hopping between the 
first and the last site of the lattice is possible. 

The parameters of the Bose-Hubbard Hamiltonian are directly connected to the 
physical lattice realized by an optical standing wave The standing wave generates 
an array of microscopic potentials - the lattice sites of the Hubbard model. A large 
intensity of the standing wave, i.e., a deep lattice potential, suppresses the hopping 
between adjacent sites and confines the atoms to individual sites. At low or vanishing 
lattice amplitudes the atoms move freely through the lattice and establish long-range 
coherence. 

The dimension D of the number basis for a bosonic system is given by 
^" A^!(J-1)! ' 

which increases rapidly with the number of particles N and lattice sites J. For fixed 
filling factor Njl = 1 the basis dimension of D = 6435 results for a system of 8 sites, 
D = 92378 for a 10-site system, and D = 1352078 for 12 sites. In the number basis 
representation the Hamilton matrix is extremely sparse, since only the nearest-neighbour 
hopping term generates off-diagonal matrix elements. Due to this sparsity the lowest 
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eigenstates can be obtained efficiently using Lanczos-type algorithms [21] for system- 
sizes up to typically I = N = 12 with standard desktop computers. The eigenstates are 
expanded in the number basis, 

W) = j:Ci^U{ni,n2,---,nj}^) (3) 

a=l 

where Cj^^ denotes the expansion coefficients of the i^-th eigenstate. These coefficients 
are obtained by diagonalizing the Hamilton matrix. 

In order to treat systems of larger size, we have introduced a truncation scheme for 
the number basis [22] • In strongly correlated regimes only a few number states contribute 
to the low-lying eigenstates, since configurations with several atoms occupying a single 
site are suppressed due to the strong interaction. For the description of the ground state, 
the basis dimension can be reduced to less than a percent without a significant loss of 
quality. In contrast, the proper description of a system of weakly interacting atoms 
requires more number states in the basis. The Hamiltonian itself provides a simple a 
priori measure for the importance of number states via its diagonal matrix elements. 
Only those states with diagonal matrix elements smaller than a certain truncation energy 
EtTunc, which depends on the Hubbard parameters, are included in the basis: 

-Etrunc > {ni, - ■ ■ ,ni\B.\ni, - ■ ■ ,nj). (4) 

By adjusting the truncation energy, one can tune the basis for a more accurate 
description of the physical system or a smaller basis dimension. Detailed benchmarks 
and applications of this basis truncation scheme can be found in j22j . 



2.2. Time evolution 

Different experimental schemes have been developed to probe the dynamic response of 
atomic gases in optical lattices. One possible scheme involves the tilting of the lattice 
to induce a static force [HE]. Another method is Bragg spectroscopy [7], which has the 
advantage of not involving any side effects like Bloch oscillations or Zener tunnelling 
[3 in]. It therefore allows a very precise determination of the excitation spectrum. 

In the case of atoms in an optical lattice this method can be implemented by a 
temporal modulation of the lattice potential Via_t{x) with frequency u. A regular one- 
dimensional optical lattice generated by an optical standing wave is given by 

V,^t{x)=Vosm\kx), (5) 

with amplitude Vq and wavenumber k. The modulation of the lattice is achieved by 
introducing a time- dependent factor: 

Vi^t{x,t) = Vo [1 + Fsm{Ljt)] sm\kx). (6) 

The amplitude of the potential © is oscillating around the background value Vq with 
frequency u and amplitude VqF with F <C 1. This introduces two sidebands with 
frequencies ±u; and defines the corresponding excitation energy u [3]. 
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As mentioned in the previous section, the physical lattice (0) enters the Bose- 
Hubbard model via the parameters J, U, and Ci In order to obtain the 

time-dependent expressions of these parameters, the localized Wannier states are 
approximated by Gaussians of width a [H]. The optimal value of the width a{t) is 
determined by an energy minimization using the lattice potential © at a given time t. 
Computing the matrix elements of the kinetic energy and the interaction part of the first 
quantized Hamiltonian within this Gaussian approximation leads to the time-dependent 
Hubbard parameters: 

J{t) = Joexp[-Fsm{ujt)] , 
U{t) = f/o[l + Fsin(cut)]i/^ 

In lowest-order approximation, the temporal change of the on-site energies is directly 
given by the change of the potential Vig,t{x,t), i.e., 

ei{t) = e^fl[l + Fsm{ujt)]. (8) 

By substituting the static parameters of the Hamiltonian with the dynamical ones 
we obtain the time- dependent Bose-Hubbard Hamiltonian. 

For our simulations of the time evolution in modulated optical lattices, we start 
with an initial state given by the ground state of the Bose-Hubbard Hamiltonian for 
given parameters Jo, Uq, and eifi. Starting from this initial state, we then evolve in time 
steps At while modulating the lattice with a fixed frequency u and amplitude VqF. The 
frequency u defines the probe energy and the response of the system is measured by 
evaluating observables using the time-evolved state IV';^)- 



2.3. Numerical methods 

The time evolution is performed either by the Crank-Nicholson scheme (CN) or by a 5th 
order predictor-corrector method| (PC), depending on the system parameters. Since the 
Crank-Nicholson scheme is an implicit method which requires the solution of a set of 
linear equations§ of the basis dimension at each timestep, it is feasible only for small 
systems. Nevertheless, due to its unconditional stability and in combination with the 
basis truncation, it is a valuable tool to simulate atomic gases in the strongly correlated 
regime. 

The predictor-corrector method is an explicit scheme which allows to treat larger 
system sizes. The drawback of this method is its numerical instability, especially in 
cases where small expansion coefficients of the states are involved. Numerical errors 
in these coefficients accumulate and lead to a collapse of the computation after a few 
steps. Nevertheless, it is possible to apply this method in the weakly interacting regime, 
where only a few small coefficients appear. In connection with the basis truncation, 
which predominantly discards these number states with small coefficients, the stability 
is improved. It is, however, not feasible to evolve systems in the strongly interacting 

X Wc use the 5th order Adams-Bashforth predictors and the Adams-Moulton eorrectors. 

§ We employ the PARDISO solver [231 EH US] which comes with Intel's Math Kernel Library. 
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regime using the PC method. The initial states consist only of a few dominant number 
states, i.e., most of the coefficients in the number representation are very small. In these 
cases the truncation cannot be used to improve the situation, because one would have 
to neglect a large number of states and, consequently, would miss some features of the 
excitation spectrum. 



3. Bose gases in an optical lattice 



3.1. Setup and linear response analysis 

As a first benchmark, we study the dynamical properties of a Bose gas in a one- 
dimensional optical lattice with periodic boundary conditions. The system consists 
of = 10 bosons at / = 10 sites of a regular lattice (ej = 0). The Bose-Hubbard 
Hamiltonian of such a system is given by 



U 

Mt) = -JY1 (al+iSi + h.a.) + IT 

1=1 ^ i=i 



n, - 1] 



(9) 



In addition to the exact time evolution, we examine the excitation of the Bose 
gas using a linear approximation of this Hamilton operator as introduced in refs. 
jHl ^1 ^1 The Hamiltonian can be written in terms of the hopping and 

interaction operators, Hj and H(/, respectively, 

H = -J{Vo)Rj + U{Vo)fiu (10) 

with the time-dependent amplitude of the physical lattice potential (jH} 

Vo{t) = Vo[l + F sm{ujt)]. (11) 

We now linearize the Hamiltonian (fTUI) with respect to the perturbation by retaining 
only the lowest-order terms of a Taylor expansion in the modulation amplitude F, 



H. 



lin 



Hr 



OF 



(12) 



F=0 



This Hamiltonian depends linearly on the temporal variation of the lattice amplitude. 
The initial Hamilton operator Hq is given by Q at the time t = 0. Evaluation of the 
derivative of the Hamiltonian ()10|1 using the dependence of the Hubbard parameters on 
the lattice amplitude results in 



Hii„(t) = Ro+FVo sin(cut) 



dlnU 



dVn 



Hr 



J 



din J 



F=0 



dVa 



d\nU 



F=0 



dVa 



H, 



,(13) 



which consists of the unperturbed Hamiltonian Hq and a linear perturbation part. The 
first term of the perturbation is proportional to the unperturbed Hamiltonian and the 
second term includes the hopping operator Hj. Due to the small modulation amplitude 
the first part of the perturbation leads to a small energy shift which can be neglected. 
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Figure 1: Lowest energy eigenvalues for a system with = 10 bosons and / = 10 sites at 
interaction strengths Uq/ Jq = 30, 15, and 3, computed by diagonahzation of the Hamilton 
matrix in a truncated number basis (see text). The vertical lines mark the strongest matrix 
elements | ( | H j | ) | between ground state and excited states. The strength of the matrix 
elements is represented by the gray level of the lines, where darker lines correspond to stronger 
values. 



The second part generates excitations, since the hopping operator connects different 
eigenstates with the couphng parameter 

) • (14) 

F=0/ 

In order to identify possible excitations of the ground state |0) to higher- lying states 
\v) of the static spectrum, we look for non- vanishing matrix elements |( z/ | Hj |0)| of 
the hopping operator Hj. 

As a first application and a preparation for the discussion of superlattice potentials 
we examine a system of 10 bosons in a regular lattice with 10 sites. The low- lying 
eigenstates are obtained using a basis consisting of the 5911 energetically lowest number 
states out of the complete basis {D = 92378), generated for the interaction strength 
Uq/Jq = 30 with the truncation energy i?trunc/<^ = 90. Figure [T] shows the spectra of 
the truncated system for interaction strengths Uq/Jq = 30, 15, and 3. 

Deep within the Mott-insulator regime, for Uq/Jq = 30, we obtain the characteristic 
gapped energy spectrum shown in figure ^a). The vertical lines illustrate the strongest 
hopping matrix elements | ( | H j | ) | between ground and excited states. The energy 
scale is shifted with respect to the ground state energy, so that the vertical lines represent 
excitations into the Hubbard bands (t/o-band) and the corresponding energies. It is 
remarkable that there are no sizable transition matrix elements to states in the 2Uq- 
band but to the 3t/o-band, as was pointed out in reference [HI [121 • With decreasing 
interaction strength the gaps are reduced and eventually vanish in the superfluid regime, 
as the sequence of plots in figure [T] demonstrates. 

A simple interpretation of the excitations to the individual Hubbard bands can 
be given in the language of particle-hole excitations. In the limit of strong repulsive 
interactions, i.e., deep in the Mott-insulator regime, the ground state of a system 



K = — JFVo sin(ci;t) 



dlnJ 



dlnU 
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schematic 




type 


energy transfer 


n A n ^ 


n. n 


one-particle-one-hole (Iplh) 


Uo 




n. £L 


two-particle-two-hole (2p2h) 


2f/o 


Ox 

n ± n 1. 


oo| 


three-particle-three-hole (3p3h) 


3f/o 


n. n 


2p2h with two particles at same site 


3f/o 



Table 1: Basic types of particle-hole excitations of first order, sorted by energy transfer. 



with unit filling is given by |0) ~ |1, 1, 1, 1, 1, 1, . . .). The simplest mechanism to 
excite this state is a one-particle-one-hole excitation (Iplh) as illustrated graphically 
in table [H The excitation energy associated with such a process equals the interaction 
strength AE = Uq, hence it represents a possible excitation from the ground state to a 
state in the first Hubbard band. Table^hsfs the basic particle-hole excitations together 
with the Hubbard band they connect to. 

3.2. Time evolution 

In order to probe the dynamical behaviour of the Bose gas we choose a fixed ratio of the 
interaction to tunnelling strength Uq/ Jq. The ground state of the system is obtained by 
exact diagonalization and is used as the initial state for the time evolution. Starting from 
this state the system is evolved in time while the lattice is modulated with frequency 
uj and amplitude FVq = 0.1 Vq. The response of the system is probed via the energy 
transfer evaluated at each timestep using the time-evolved state |^E',t): 

AE(t) = (^,t|Ho|*,t)-^o. (15) 

Ho is the Hamiltonian at time t = 0, and Eq is the ground state energy of the initial 
Hamiltonian. 

Figure |2l shows the energy transfer for 10 bosons and 10 lattice sites at several 
ratios Uq/Jq, ranging from the Mott insulating phase (MI), depicted in panels a)- (e), 
to the superfiuid regime (SF), shown in panel |21^f). The density plots in the lower part 
of each panel show the energy transfer AE as a function of the modulation frequency uj 
and time. The line plots in the upper part of each panel show the energy transfer as a 
function of u and averaged over the full time evolution. The (red) arrows indicate the 
energies of excited states with sizable hopping matrix elements to the ground state, the 
size of the matrix elements is reflected by the length of the arrows. 

In the simple picture of particle-hole excitations (see table one would expect 
excitations at integer multiples of the interaction strength in the strongly correlated 
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Figure 2: Energy transfer /S.E for a system of = 10 bosons and / = 10 lattice sites at 
different interaction strengths due to the modulation of the lattice amphtude {F = 0.1). The 
density plots show the energy transfer as function of modulation frequency and time, the 
line plots on top depict the energy transfer as function of frequency averaged over the whole 
evolution. The sequence from (a) to (f) corresponds to interaction strength Uq/Jq = 30, 20, 
15, 10, 5, and 3. The (red) arrows mark the energies of eigenstates {v) which are connected 
to the ground state by the hopping operator Hj (compare figure 
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Figure 3: Relation between the excitation energies of the static Hubbard Hamiltonian and the 
resonance structure of a bosonic system (I = N = 10) at interaction strength Uq/Jq = 20. 
The left panel shows the time-averaged energy transfer as a function of the frequency w; the 
right panel depicts the lower part of the eigenvalue spectrum of the corresponding static 
Hamiltonian. The vertical gray lines mark the eigenstates which are connected to the 
ground state by strong matrix elements |(i^|Hj|0)|. The horizontal (red) lines point out 
the correspondence between the excitation energies and the resonance structure. 



regime. The simulations illustrated in figure |21I^a)-(d) confirm this assumption partially. 
A large resonance appears at the frequency u = Uq in the Mott phase, in panels (b)-(d) 
weaker excitations also appear at u; = SUq. In the case of the simulation for Uq/ Jq = 30 
the truncated basis does not include number states with energies of 2Uq and higher, and 
consequently the 3[/o-resonance does not appear. 

The absence of a resonance peak at frequency u = 2Uq was already evident from the 
linear response analysis illustrated in figure ^a)-(c), which shows no significant matrix 
elements between the ground state and excited states in the 2f/o-band. This reveals that 
2p2h-processes as schematically depicted in table ^ are not relevant for the excitation. 

Towards the superfluid regime the resonances are significantly broadened, which 
is in agreement with experimental evidence [H] . Figure Efe) shows a simulation in the 
transition region from Mott-insulator to superfluid around {Uq/Jq)c = 4.65 0. The 
response spectrum reveals a significant shift of the resonance towards higher frequencies 
for decreasing interaction strength. The simple approximation of a single number 
state as the ground state is not adequate for weaker interactions and especially in the 
superfluid regime due to the delocalization of the atoms. The admixture of other number 
states leads to additional particle-hole excitations and therefore to a broadening of the 
resonances. 

The fine-structure of the resonances (figure can be explained by the linear 
response analysis. Figure El reveals the connection between the excitations to higher 
eigenstates and the structure of the resonances appearing in the time evolution. The left 




Figure 4: Comparison of the time-averaged energy transfer for a Bose gas with I = N = 8 
at the interaction strength Uq/Jq = 30. The different hnes show the result for the fuh basis 

( ) and for truncated bases with -E'truncZ-^o = 90 ( ), Etmnc/Jo = 60 ( -), and 

-E'trunc/'^o = 30 ( ). Depicted are the main resonance at uj = Uq (a), the non-Unear 

resonance at uj = Uq/2 (b), and the 3C/o-resonance (c). 



hand side shows the time-averaged energy transfer ofa/ = A^ = 10 bosonic system, the 
right hand side the corresponding energy spectrum obtained by solving the eigenvalue 
problem of the initial Hamiltonian. The vertical lines indicate the eigenstates that are 
connected to the ground state by large matrix elements of the hopping operator Hj. The 
horizontal dashed (red) lines and the (red) arrows mark the corresponding excitation 
energies, which resemble — with remarkable precision — the fine-structure of the 
resonances emerging in the time evolution. 

In addition to the resonances at integer multiples of the interaction strength, there 
is also a tiny resonance at Uq/2 in the strongly interacting cases. It results from non- 
linear effects, such as the absorption of two photons of the energy ?7o/2, which are not 
captured by the linear response analysis. 

3.3. Benchmark for the truncation scheme 

In section ini we introduced the basis truncation scheme in order to reduce the numerical 
effort of the simulations. For static ground state properties, we have shown that the basis 
dimension can be reduced typically by two orders of magnitude without a significant loss 
of precision [22] • This does not necessarily hold for the exact time evolution discussed 
in this paper. In order to describe the excitation effects properly, one has to use a basis 
which allows to describe more than a few eigenstates. 

To test the basis truncation in the context of time-dependent calculations with 
a modulated lattice we compare time evolutions with several cut-off energies -Etrunc- 
We perform the time evolution of a system of = 8 boson in J = 8 sites for 
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-E'trunc/ ^0 


basis dimension 


excitations included 


OO 


6435 (complete) 


all 


90 


1205 


up to 3-particle-3-hole 


60 


477 


up to 2-particle-2-hole 


30 


57 


up to 1-particle-l-hole 



Table 2: Main characteristics of the bases used for the comparison. Shown are the truncation 
energy -Etruno the basis dimension and the highest included particle- hole excitation with 
respect to the reference state |1, 1, 1, 1, 1, 1, 1, 1). 

Uq/Jo = 30, which is the largest system with integer filling that allows the use of 
the complete basis. We compare the energy transfer as a function of the modulation 
frequency u, averaged over the evolution time tmax/ Jq^ = 10. Table El specifies the 
characteristics of the bases used for the comparison. Figure EJa) shows the results for 
the [/o-resonance. The response obtained with the basis truncated at -Etrunc/ Jo = 90 
is in excellent agreement with the calculation using the complete basis. The basis 
with cut-off energy -Etrunc/ Jo = 60 shows small deviations in the peaks. The basis 
truncated at -Etrunc/ -^o = 30 reproduces the peak as well as its fine-structure but the 
whole resonance is shifted by roughly AE = 1.5Jo. The shift towards higher energies 
for more severe truncations can be explained by the Hylleraas-Undheim-MacDonald 
theorem [211 123 ? which states that the exact eigenenergies are lower bounds for the 
corresponding eigenvalues obtained with a truncated basis. In a variational picture, the 
basis truncation reduces the flexibility of the number state representation (trial states) 
and thus leads to an increase of the energy eigenvalues in the truncated space. 

Figure l^b) shows an enlarged energy interval around the non-linear resonance 
at frequency u = Uq/2. The basis with the truncation energy -Etrunc/^^o = 90 
perfectly reproduces the result of the complete basis, even the -Etrunc/ Jo = 60-basis 
causes only small deviations. In agreement with the results for the f/o-resonance, the 
whole resonance is significantly shifted to higher energies in the simulation using the 
-Etrunc/ -'o = 30-basis due to the overestimation of the eigenenergies. A closer look at 
the 3[/o-resonance in figure |3fc) illustrates the limitations of truncated bases for the 
description of high-lying resonances. For truncation energies Etmnc/Jo = 30 and 60 
the bases do not include any number states that correspond to energies as high as the 
excitation energy, hence no energy transfer is possible. The basis with the truncation 
energy -Etrunc/^o = 90 just about includes the numbers states matching the excitation 
energy of = SUq = 90 Jq, and consequently, an energy shift occurs, which is comparable 
in size to the shift of the t/o-resonance for the truncated basis with Etmnc/ Jo = 30. 



Response of Base gases in time- dependent optical superlattices 



13 




Uo/Jo 

Figure 5: Left panels: Two different sets of on-site energies Ci of period-five two-colour 
superlattices used in this paper. Right panel: Phase diagram of bosons (N = 10) in an 
optical superlattice (/ = 10) spanned by the interaction strength Uq/Jq and the superlattice 
amplitude A/Jq for the set (A) of on-site energies. Depicted is the maximum coefficient C^^^ 
of the number state expansion ^ of the ground state. 

4. Bose gases in a two-colour superlattice 
4.I. Two-colour superlattice and the phase diagram 

In this section the dynamical signatures of Bose gases in two-colour superlattices are 
investigated. In experiment, superlattices can be generated by a superposition of two 
optical standing waves, which leads to a sinusoidal modulation of the depth of the 
individual lattice wells. The spatial modulation enters the Hubbard model via the 
external potential term discussed in section 12.11 with an appropriate distribution of the 
on-site energies e^. The left hand side of figure shows two different distributions of 
on-site energies for a two-cell superlattice of 10 sites. Essentially, the two superlattices 
differ in the relative phase of the standing waves. We will use these two realizations to 
estimate the dependence of the response on the detailed topology of the superlattice. In 
the following we refer to these on-site energy distributions as superlattices (A) and (B). 

The amplitude A of the superlattice is an additional parameter which generates 
a rich structure in the phase diagram IT^ . The right hand side of figure El shows 
the two-dimensional phase diagram for = 10 bosons and J = 10 sites with on-site 
energies according to superlattice (A) (figure 13)) . Plotted is the maximum coefficient 
Cmax of the expansion Q of the ground state in the number basis. The darker (red) 
shadings represent regions of a large maximum coefficient which indicates that a single 
number state is dominating the ground state. Brighter shadings refer to small maximum 
coefficients which correspond to a ground state given by a superposition of many number 
states. 

The bright region at small interaction strengths and superlattice amplitudes 
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Figure 6: Schematic diagram of a single cell of the superlattice: If the superlattice amplitude 
A is close to the interaction strength Uq, the number states (a) and (b) have roughly the same 
energy. The relative height between the sites represents the different on-site energies e^. 

indicates the superfiuid phase. By increasing the interaction strength the system enters 
the homogenous Mott-insulator phase, which is represented by the dark (red) area below 
the diagonal {Uq = A). As long as the interaction strength exceeds the superlattice 
amplitude it is energetically unfavourable to occupy a site with more than one boson. 
Close to the point where amplitude and interaction are equal, the maximum coefficient 
Cmax decreases: Here, the two number states depicted in figure El have the same energy, 
since the energy gain in on-site energy obtained by hopping to the deepest well (figure 
infb)) is compensated by the increase in interaction energy. 

These redistributions of the particles in favour of the deep lattice wells become more 
and more important if the superlattice amplitude increases further. As soon as the gain 
in on-site energy for a particular redistribution exceeds the loss due to the increased 
interaction energy, the particle will move to a deeper lattice well. These successive 
redistributions generate the lobe structure in the phase diagram (figure E}- In the case 
of random lattices the redistributions happen almost continuously, giving rise to the 
Bose glass phase. Eventually all atoms are localized at the deepest lattice well. This 
complete localization also appears for small or vanishing interaction strengths Uq/Jq in 
the phase diagram. 

4.2. Linear response analysis 

As a first analysis of the response of superlattice systems we examine the spectrum of 
the static Hamiltonian (^J) for a single superlattice cell. Figure [7| shows the spectra of 
the Bose-Hubbard Hamiltonian for = 5 bosons and 1 = 5 sites at fixed interaction 
strength Uq/ Jq = 30. Instead of two-cell superlattices such as (A) and (B) we consider a 
single cell only, because the small basis dimensions allow us to solve the full eigenproblem 
without truncation. Note that the energies plotted in figure 13 are shifted by the ground 
state energy in such a manner that the ground state energy is zero. 

The sequence of plots in figureElreveals the effect of the superlattice amplitude A on 
the Hubbard band structure. The vertical lines mark the eigenstates which are connected 
to the ground state via nonvanishing matrix elements of the hopping operator Hj. The 
relative strength of these matrix elements | ( z/ | H j | ) | is indicated by the gray level of 
the lines. Darker shadings correspond to larger matrix elements. 

In the absence of a spatial modulation the typical Hubbard bands of the strongly 
correlated regime are visible, as depicted in figure [7|^a). As discussed in sectional the 
second Hubbard band (2f/o-band) is not connected to the ground state in first order. 
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Figure 7: Energy spectrum of an / = = 5 system at fixed interaction strength Uq/ Jq = 30 
and various superlattice amplitudes A/Jq. The vertical lines mark the eigenstates that are 
connected to the ground state via the hopping operator Hj. 



Figure Efb) shows the eigenenergy spectrum for a system with a small superlattice 
amplitude of A/Jq = 5. The gapped structure is still visible but the width of the 
individual bands is significantly increased and the gaps are reduced. The breaking of 
the spatial symmetry by the superlattice generates a larger number of matrix elements 
which connect the ground state to higher bands. At a superlattice amplitude A/ Jo = 10, 
shown in figure Hl^c), the gaps have completely vanished, but nevertheless, the largest 
matrix elements still couple the ground state to eigenstates in the Uq and 3t/o bands. 

At an amplitude A/Jq = 20, which is still within the homogenous Mott-insulator 
phase, the separation of the sizable matrix elements which were initially associated with 
the Uq and 3?7o bands dissolves as shown in figure [Z|^d). An interval of large matrix 
elements remains for excitation energies up to 1.5t/o; which resembles the original Uq- 
band. Beyond that, there is a number of weaker matrix elements covering excitation 
energies of up to 3.5f/o. With increasing superlattice amplitude, any semblance of 
the band structure disappears. As depicted in figure [7|^e) for A/Jq = 50 one finds 
a continuous distribution of matrix elements with dominant matrix elements at low 
excitation energies. For very large superlattice amplitudes, e.g. for A ^ 2Uq as 
depicted in figure d^f), the higher- lying matrix elements are suppressed and strong 
matrix elements appear only at the lower end of the energy scale. This characteristic 
behaviour emerges also from the dynamical simulations. 
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Figure 8: Energy transfer during lattice modulation (F = 0.1) for a superlattice system with 
/ = = 10 and interaction strength C/q/Jq = 30 for two superlattice amplitudes A < C/q in 
the homogenous Mott-insulator phase. The line plots show the time-averaged energy transfer 

as function of m for the superlattices (A) ( ) and (B) ( ) defined in figure 101 The 

(red) arrows mark the excitation energies predicted by linear response analysis. The density 
plots illustrate the energy transfer as a function of frequency and time for the superlattice (A). 



4-3. Time evolution 

We now turn to fully dynamical simulations of Bose gases in time-dependent 
superlattices. We focus on a system of = 10 bosons in the two-cell superlattices with 
/ = 10 sites as defined by (A) and (B) in figure AH simulations are performed for a 
fixed interaction strength Uq/Jq = 30 and various values for the superlattice amplitude 
A/Jq. As initial state we use the ground state of the static Bose- Hubbard Hamiltonian 
for the same parameters. In analogy to the procedure for regular potentials in section 
13 the lattice potential is modulated in time with a frequency u and a fixed relative 
amplitude F = 0.1. 

Figures |H1 and IHl show the energy transfer to the system at several superlattice 
amplitudes for a fixed interaction of Uq/Jq = 30. All simulations are performed with 
each of the two on-site energy distributions depicted on the left of figure El in order to 
assess the impact of the detailed distribution of the on-site energies e^. The density 
plots in the lower part of each panel show the energy transfer as function of time and 
frequency for superlattice (A). The plots in the upper part represent the time-averaged 
energy transfer as function of the frequency for both superlattice topologies. The (red) 
arrows above the individual peaks mark the excitation energies associated with the 
strongest matrixelements |(z/|Hj|0)|. The size of the arrows indicates the relative 
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Figure 9: Energy transfer during lattice modulation (F = 0.1) for a superlattice system with 
I = N = 10 and interaction strength Uq/Jq = 30 for three superlattice amplitudes A > C/q in 
the quasi Bose-glass phase. The line plots show the time-averaged energy transfer as function 

of uj for the superlattices (A) ( ) and (B) ( ) defined in figure The (red) arrows 

mark the excitation energies predicted by linear response analysis. The density plots illustrate 
the energy transfer as a function of frequency and time for the superlattice (A). 

strength of the matrixelements. 

Figure |H^a) depicts the t/o-resonance in the case of a regular lattice, i.e., for 
vanishing superlattice amplitude A/Jq. As pointed out before, the excitation energies 
associated with the strongest matrix elements resulting from the linear response analysis 
nicely describe the position and fine-structure of the resonance. If we increase the 
superlattice amplitude to the value A/Jq = 20 — still remaining in the homogeneous 
Mott-insulator phase — the overall width of the resonance structure increases. The 
characteristic scale of the fine structure increases as well, and additional peaks emerge. 
This is in accord with the results of the linear response analysis of the previous section. 
An increase of the superlattice amplitude leads to a broadening of the Hubbard bands 
(cf. figure [7| (b)) which corresponds to the broadening of the resonance. The larger 
number of peaks in the resonance corresponds to the larger number of matrix elements. 
The comparison of the two different superlattice topologies shows that the envelopes 
of the resonance are practically identical in both cases. Only the details of the fine 
structure depend on the particular set of on-site energies used. 
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After crossing the transition to the quasi Bose-glass phase at A = [/q, the resonance 
structure changes dramatically as shown in figure IHl Already at A/ Jq = 35, i.e., slightly 
above the transition at A/Jq = Uq/Jq = 30, the resonance is completely fragmented. 
There is no longer a smooth envelope with centroid at u = Uq. Instead, the strength 
is split into two groups of individual peaks: one group around the original resonance 
position and another group at low excitation energies. In particular, we observe low- 
lying resonances atu;/Jo < 10 — a regime where no response was observed for a system 
in the homogeneous Mott-insulator phase. With increasing superlattice amplitude, the 
response is shifted towards lower excitation energies, as indicated by the sequence of 
plots in figureini At A/Jq = 50, for instance, no response is left at the original resonance 
position u = Uq and all peaks are concentrated at low excitation energies. 

This characteristic behaviour of the response appears to be a clear signature for the 
Mott-insulator to quasi Bose-glass transition for boson in superlattices, which is directly 
accessible to experiments. The comparison of the two different superlattice topologies, 
both with a period of five sites, demonstrates that the fine-structure of the response 
depends on the details of the superlattice, but that the gross characteristics are not 
affected. 

Finally, one should note that the linear response analysis already hints at these 
substantial changes in the resonance spectrum. For larger superlattice amplitudes 
however, effects beyond the simple linear perturbation scheme become increasingly 
important and lead to significant discrepancies in comparison to the full time-dependent 
simulation (cf. figure IHfc)). 

5. Summary and conclusions 

We have studied the dynamics of Rose gases in one dimensional optical lattices and 
superlattices with time dependent lattice amplitudes. The time evolution based on the 
time-dependent Bose-Hubbard Hamiltonian is performed numerically for systems with 
up to 10 bosons and 10 lattice sites. In order to reduce the numerical effort we have used 
an a priori basis truncation scheme [22], which reduces the dimension of the number 
basis to a tractable size. The control parameter of this truncation scheme is the cut-off 
energy -Etrunc? which defines the maximum energy of the number states in the basis. We 
have shown that the truncation also allows reliable dynamical calculations. 

As a first application, we have examined Rose gases in a regular lattice potential. 
In agreement with experiment jH] and other theoretical results jHl EI E] , we illustrated 
the characteristic resonance structure of the Bose gas in the strongly correlated Mott- 
insulator regime, which is washed out and broadened towards the superfiuid phase. 

Treating the time-dependence of the Bose-Hubbard Hamiltonian as a linear 
perturbation demonstrates that excitations are generated by the hopping operator in 
first order pT H ITT | IT^. The analysis of matrix elements of the hopping operator between 
ground and excited states allows a prediction of the resonance position and a detailed 
explanation of the fine structure. It is shown, that the individual peaks of a resonance 
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correspond to large matrix elements |(0 | Hj | where {u) are the eigenstates of the 
initial Hamiltonian. 

In the second part we have investigated dynamical signatures of quantum phase 
transitions in two-colour superlattice potentials. We have shown that for superlattice 
amplitudes A smaller than the interaction strength Uo, i.e., in the homogeneous Mott 
insulator phase, the resonances resemble those of regular lattices. With increasing 
superlattice amplitude, the resonance structures become broader and additional peaks 
appear in their fine-structure. This finding is consistent with recent experiments 
[T7] using time-dependent superlattices with incommensurate wavelengths of the 
superimposed standing waves. First calculations for incommensurate superlattices 
confirm that the gross characteristics of the resonance spectrum are independent of 
the precise lattice topology, only the fine-structure of the resonances is affected. 

As soon as the superlattice amplitude exceeds the interaction strength, i.e., on 
entering the quasi Bose-glass phase, the resonance structure changes completely. The 
resonance a.tu = Uo is fragmented and a strong low-energy component develops. Further 
increase of the superlattice amplitude eventually suppresses all strengths at the original 
resonance position, such that only the response at very low excitation energies remains. 
This characteristic behaviour might serve as an experimental indicator for the transition 
from homogeneous Mott-insulator to quasi Bose-glass phase. 

In order to address the experimental scenario ^\ directly, we are going to 
investigate the influence of different lattice topologies, external trapping potentials, and 
filling factors in detail |2H] • Nevertheless, the results presented in this paper apply to the 
Mott domain with one atom per site even in the presence of a weak trapping potential. 
Another topic for future studies is the impact of non-zero temperatures on the response 
behaviour, which is beyond the present zero-temperature formalism. 
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